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We investigate 4-quark {qqqq) systems as well as multi-quark states with a large number of 
quarks and anti-quarks using the chromodielectric model. In the former type of systems the flux 
distribution and the corresponding energy of such systems for planar and non-planar geometries are 
studied. From the comparison to the case of two independent gg-strings we deduce the interaction 
potential between two strings. We find an attraction between strings and a characteristic string flip 
if there are two degenerate string combinations between the four particles. The interaction shows 
no strong Van-der-Waals forces and the long range behavior of the potential is well described by a 
Yukawa potential, which might be conflrmed in future lattice calculations. The multi-quark states 
develop an inhomogeneous porous structure even for particle densities large compared to nuclear 
matter constituent quark densities. We present first results of the dependence of the system on the 
particle density pointing towards a percolation type of transition from a hadronic matter phase to 
a quark matter phase. The critical energy density is found at Sc = 1.2 GeV/fm^. 
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I. INTRODUCTION 

Quantum Chromodynamics (QCD) is the theory for 
color charged quarks and gluons. This theory has been 
tested successfully in the regime of large momentum 
transfer such as in deep inelastic scattering, where per- 
turbative methods can be used due to asymptotic free- 
dom. It is a common belief that QCD should be able to 
describe all systems ruled by strong interactions. These 
cover a wealth of different regimes ranging from the dy- 
namics of quasi free quarks and gluons in a quark gluon 
plasma (QGP) at high temperatures or densities, over 
the formation of hadrons out of quarks to the interac- 
tions between those color neutral hadrons. However, in 
this latter region of small relative momenta, where the 
confinement phenomenon plays a dominant role, no con- 
vincing analytical techniques have been established yet. 
The only calculations from first principles are restricted 
to QCD lattice simulations, where the results are limited 
due to today's computing power. 

Therefore one still has to rely on models, that de- 
scribe the interactions between quarks and gluons bound 
to hadrons phenomenologically and that include confine- 
ment. Such models are for example the well known MIT 
bag model P, Q , where quarks are confined by a given 
external cavity or the quark molecular dynamics model, 
[3-5], where quarks follow the Hamiltonian dynamics 
subject to a linear rising potential between quarks (q) 
and anti-quarks (q). Another approach is the model of 
the dual super conductor known also dual Abelian Higgs 
model or dual Ginzburg-Landau model ^ iTj where con- 
finement is achieved by monopole condensation [^|^ and 
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an accompanying magnetic supercurrent. The model of 
the stochastic vacuum relies on the calculation of Wilson 
loops in a Gaussian approximation [lOl - fT^ which leads 
to a linearly rising gg-potential. 

In this work we choose the framework of the chro- 
modielectric model (CDM) [l^jl^. Unlike the MIT bag 

model, CDM has the benefit, that bags with a smooth 
surface are created self-consistently and dynamically out 
of the underlying field equations, due to the presence of 
colored quarks. 

The CDM has already been used to calculate prop- 
erties of the nucleon and its low lying resonances like 
masses, magnetic moments and the axial- vector/vector 
coupling constant ratio jl^]. In (ITrilOj a description of 
qq-strings was given, and in poj the parameters of the 
CDM were adjusted to reproduce results of lattice calcu- 
lations J.1, 22] . In the same work, the fiux tube struc- 
ture of a baryon like qqq-hag has been studied. Further, 
within the model, the interactions between gq-strings 
and the nucleon-nucleon interaction in vacuum and in nu- 
clear matter has been discussed. In another ap- 
proach the model has been used in a transport theoretical 
framework to describe the dynamics of quarks bound in 
nucleons and strings [2^l30l ]. In "sT] the disintegration of 
qq-pairs in strong color electric fields has been observed. 
A full molecular dynamical simulation of hadronization 
out of a gas of quarks and gluons has been presented in 

M 

In this paper we will analyze the interactions between 
the color electric flux tubes for a wide class of different 
quark configurations. In a previous paper [20] we have 
given the structure of meson like qq-states and baryon 
like qqq states. The model has been successfully adjusted 
to reproduce lattice results for both the qq-potential and 
the transverse shape of the flux tubes. It is an open 
question, how the flux tubes of the basic white two- and 
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three-particle clusters interact with each other. In order 
to understand the interactions we extend our previous 
analysis to configurations with more than three particles. 
The easiest of such systems, the qqqq-system, already 
develops two distinct bags, that may interact with each 
other. We can study those systems for different spatial 
quark configurations. Besides this, the system is still sim- 
ple enough to be treated on the lattice. In fact there are 
lattice results for the four-particle system in SU(2) 
and also in SU(3) [s^- Where possible, we will compare 
our results with those obtained on the lattice. Knowing 
the interactions between the flux tubes of color neutral 
objects, it is an interesting issue, how this interaction 
governs the transition from a system of distinct white 
hadrons to a system of interacting colored quarks in a 
quark plasma. It is an old prediction of QCD, that there 
is a rapid transition for increasing temperature |35| . This 
transition to the quark gluon plasma (QGP) has been ex- 
plicitly seen for vanishing baryon chemical potential as a 
steep rise of the thermodynamic pressure of the system 
at a temperature T = (170-180) MeV For non-zero 
baryo-chemical potential, lattice calculations suffer from 
technical problems and cannot be easily performed up 
to now. See [s^ l38|. where the transition temperature 
was explored for small chemical potential. However, for 
increasing baryon densities one expects that the hadrons 
start to overlap and the quarks are free over a much larger 
volume and a transition to a quark gas might occur as 
well. In contrast, in our treatment of the CDM we have 
only static configurations and therefore we do not de- 
scribe quark systems at non-zero temperature. But it is 
possible to vary the quark density and study many-quark 
systems and the behavior of the very dense flux tubes. 

The paper is organized in the following way. In sec- 
tion^ we introduce briefly the chromodielectric model, 
give the equations of motion for the underlying fields as 
well as the corresponding field energies. We solve these 
equations numerically in three spatial dimensions. We 
also discuss the color structure of the model and its con- 
nection to the SU(3) color algebra. Section imi is devoted 
to the long range behavior of fields of gg-strings, which 
extends the discussion of the bulk properties of strings 
in [23. We will show the characteristic exponential de- 
cay of all fields with the distance from the string center. 
In section HVI we show the interactions between two qq- 
strings of various lengths and with relative orientations 
to each other. We first show our results for the color 
electric fields in section HV Al and the distribution of the 
energy density in section HVBI The distortion of the one 
qq-string under the influence of the other one is shown ex- 
plicitly. For simplicity, we concentrate on the discussion 
of strings that are parallel or anti-parallel to each other. 
In this case, all particles are lying in a single plane. How- 
ever, as the calculations are done in three dimensions, 
we are able to study particle configurations too, that are 
extended in three space dimensions. We will show the 
results of the calculations for these tetrahedronic config- 
urations as well. In section llV CI we give the interaction 



potential of the two strings as calculated from field en- 
ergies. Again we concentrate on the potential for plane 
configurations but we will discuss also a wide range of 
3-dimcnsional configurations. We will extract the static 
string interaction potential which can be compared to 
lattice calculations. In section [3 we present first results 
on static multi-particle systems. We show the heteroge- 
neous structure of the emerging flux tubes and analyze 
the scaling of the multi-quark properties with the particle 
density. A qualitative description of the deconfinement 
phase transition is given. Finally we discuss our results 
in section IVll and give our summary. 



II. THE CHROMODIELECTRIC MODEL 

As there is no confinement in quantum electrodynam- 
ics, one believes that confinement is due to the non- 
Abelian nature of QCD. Although the detailed mecha- 
nism of confinement has not been revealed within QCD, 
there is strong evidence for it from lattice QCD. The most 
prominent result is the linear rising potential between a 
static quark anti-quark pair at zero temperature [3^ l40| 
and at finite temperature In addition, the forma- 

tion of long, string like flux tubes between a qq-pair has 
been seen both in SU(3) [Ulil and in SU(2) where 
the results are much more precise for numerical reasons. 
Despite the fact that the structure of the QCD vacuum 
is highly complicated, its long range behavior might be 
transparent. In the CDM it is assumed that this vac- 
uum behaves as a perfect dielectric medium described by 
a vanishing dielectric constant Kvac- Dual to a normal 
superconductor, where magnetic fields are expelled from 
the superconducting phase, in CDM the color electric 
fields are expelled from the QCD vacuum. In the pres- 
ence of color charged quarks, the resulting color fields are 
compressed into flux tube like excitations of the vacuum. 
In [24l l44l I45I I a renormalization group derivation was 
given for the lattice colordielectric model, which has a 
scalar confinement field, and keeps strongly coupled non- 
Abelian fields in the large distance limit. CDM is formu- 
lated in terms of two Abelian color fields A'^''^ only and 
an additional scalar confinement field a. This scalar field 
is designed to already include the non- Abelian effects of 
the gluon sector. There are indications that in the long 
range limit only the Abelian components contribute to 
the observable quantities: 't Hooft suggested in Q the 
maximal Abelian gauge for projecting out a Cartan sub- 
group believed to be relevant for infrared aspects of QCD. 
In 46.' 48j further support for the Abelian dominance was 
found due to the mass generation of off-diagonal gluons. 
In |49l Isoj the string tension in the Abelian approxima- 
tion was reproduced within some percent deviation from 
the full value. 

Following this reasoning, CDM is formulated in a Car- 
tan subgroup of QCD reducing the independent color 
fields to a set of two commuting field A^^" with a e {3, 8}. 
These two fields are connected to the Gell-Mann matri- 



ces A'^'^ which commute with each other in the standard 
representation. Confining effects are merged into the di- 
electric coupling of these gluon fields to the dielectric 
medium generated by the confinement field a. The CDM 
Lagrange density can now be given as 



C = 



1^9 



a e {3,8} 



(la) 
(lb) 
(Ic) 
(Id) 



where A^-"" — (cj)"', A"") is composed out of the color elec- 
trostatic potential cj)'^ and the vector potential 

In the following we are interested especially in the in- 
teractions of the electric flux tubes, thus we have omit- 
ted the dynamic term for the quark degrees of freedom. 
Quarks enter into the model only via the external color 
current — (p", J") with a coupling strength gg. Fur- 
thermore we treat the quarks as infinitely heavy, static 
sources. Therefore the color current vanishes J*" = 0. 
The color charge density p°(a;) = J^k 9k'^i^~^k) is given 
as a sum over all quarks with two color charges g° and 
a spatial distribution w{x). In principle the quarks are 
point-like objects, but for numerical reasons we assign a 
Gaussian distribution w{x) = (27rrQ)~^/^ exp(— af^/2rQ) 
with a small width vq — 0.02 fm, which is on the order 
of the spacing of the numerical grid used in the calcula- 
tions. The grid spacing chosen is small compared to all 
physical sizes like the flux tube radius. 

In the Abelian approximation the quarks still have 
three different colors. These are expressed as three di- 
mensional unit vectors in the fundamental representation 
of color space as \c) £ {\r) , \g) , \b)} . The color charges are 
then given by the diagonal entries of the corresponding 
generators t°- — A"/2, i.e. = (c|t"|c). Formulated dif- 
ferently, the two color charges {q'^,q^) are given by the 
weight vectors of QCD jj. The numerical values of the 
color charges can be read off from tableland are depicted 
in figure n 
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TABLE I: The color charges of the three colors with respect 
to the two Abelian color fields. The charge of the anti-colors 
are given by the negative of the color charges. 



The model inherits a U(l)xU(l) gauge symmetry and 
a global symmetry corresponding to finite rotations in 
color space |23| which is a remnant of the original SU(3) 
symmetry of QCD. Only the Lagrange density and the 
energy density, which is stated explicitly later on, are 
invariant under these rotations but not the color fields 
themselves. 
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FIG. 1: The color charge (q^, 2\/3g*) with respect to the color 
fields A^'^ for the three colors r,g,b. 



The color field tensor F^"'^'' (a G {3, 8}) in Eq. ||Tdl de- 
fines the color electric and magnetic fields E°- ~ _^Oi,a _ 

(-V(/)'' - 9*1°) and Bf = -^le^jkF^^^" = (V x A% 
respectively. With the definition of the medium fields 
= K{a)E°- and i?" — K{a)B°' the equations of motion 
following from Eq. are given by 

V-^" - gsp^ , (2a) 

V X i?-^ - dtD- = " , (2b) 

V X + dtB" = , (2c) 
V-B = , (2d) 

a + U\a) = -\K'{a)F;^F^-'^ . (2e) 



With J° = and the assumptions that all time deriva- 
tives vanish exactly and A° = 0, the magnetic field is 
B"^ = and the equations of motion (c.o.m.) can be cast 
into the following form: 



V • (K(cr)V<?!>") = 

VV-C/'(ct) = 



-9sP , 
1 K'{a) 
'2^2 (a) 



(3a) 
(3b) 



The energy of the system for static configurations is given 
by 



-E'tot — Eel + Eyoi + Esu, 
d^rE" -D" 



Ef,] — 



E, 



vol 



d?rU{(T) 



(4a) 
(4b) 

(4c) 

(4d) 
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where we have labeled the different energy parts as to- 
tal energy, electric energy, volume energy, and surface 
energy, respectively, as explained in 0. 

The confinement field is exposed to a quartic self in- 
teraction 



U{a) = B + aa^ + ba^ + ca* 



(5) 



Two specific forms with different parameters B, a, b, c are 
shown in fig. |2] The generic form of U develops two 
(quasi-) stable points, which separates the two distinct 
phases of the model. The first defines the vacuum ex- 
pectation value of the scalar field a = dvac and is asso- 
ciated to the energy density of the confined phase with 
U{a) = 0. The second at tr = is associated to the 
deconfined phase with an energy density U{0) ~ B. We 
refer to the former phase as the non-perturbative vacuum 
and to the latter as the perturbative vacuum. Given the 
generic form of the potential U and the vacuum value 
(Tvac , there are only two other independent parameters de- 
scribing U. We choose the perturbative value B = U{0) 
and the curvature nig — U" {a)\a^^^ of the potential at 
fvac, where the primes denote derivatives with respect 
to tr. nig behaves as the mass of the confinement field 
a and can be interpreted as the mass of a glueball, as 
all non-perturbative gluonic effects are collected in the 
scalar field of our model. Another possible interpreta- 
tion is to relate the mass of the confinement field to the 
mass of the off-dia gona l gluons generated in the maximal 
Abehan gauge |4al48| . The parameters a,b,c appearing 
in Eq. Q can be expressed as 



1 m^gj, - 12B 



8B 



1 ^gCTvac 



6B 



(6a) 
(6b) 
(6c) 



with the additional constraint, that a > to ensure that 
there is no relative maximum at ct = 0. 

The perturbative and the non-perturbative phases dif- 
fer not only with respect to the corresponding energy 
densities U{0) = B and C/(CTvac) = 0, respectively, but 
also in their dielectric behavior. In the former, the di- 
electric constant K{a = 0) = 1 allows for freely propagat- 
ing fields, whereas in the latter K(crvac) — ^vac ^ 1 and 
the electric fields are suppressed. Kvac = would lead to 
perfect screening of the color fields and the non-zero but 
small value of Kvac is introduced for numerical reasons 
[20|. The dielectric function k((t) is designed to interpo- 
late smoothly between the two values and we choose the 
following parameterization: 
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FIG. 2: The scalar potential U(a) 
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for the two different pa- 
I and stated explicitly in 
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with s — cr/(Tvac and with coefficients 

= 

In the limit Kvac ~^ the non-perturbative vacuum be- 
haves as a perfect dielectric medium and all electric fields 
are expelled out of regions where a — fivac- Note, how- 
ever, that Kvac = is numerically not feasible. As dis- 
cussed in [201 there are different possibilities to parame- 
terize k((t). With the polynomial form chosen here, the 
results as shown below do not depend on the exact value 
of Kvacj once it is smaller than Kvac < 10~^. In this work 
we choose Kvac = 10^^. 

The fields and a and the corresponding field en- 
ergies are calculated numerically according to Eqs. © 
and |0J| using the multi-grid FAS-algorithm given in 
aail,52j. 

The dependence of the numerical results on the model 
parameters was studied in detail in Different sets 

of parameters were given that reproduce the Cornell po- 
tential for a qq-Tpa.iT and simultaneously the profile of the 
color electric string. The Cornell potential was found in 
heavy quarkonium spectroscopy l53-|56l and was repro- 
duced on the lattice both in SU(2)0fand in SU(3) 
and is given by 



V,{R)^2CfEo-Cf- 
H 



tR 



(9) 



(7) 



We focused on reproducing the string tension r ~ 
980 MeV and the Coulomb coefficient a — 0.3. Here R is 
the distance between the quark and the anti-quark and 
Cf = 1/3 is the corresponding color factor in the Abelian 
approximation. The constant and finite term Eq is due 
to the non-zero width of the particles. The transverse 
profile was calculated on the lattice |2l|. in the frame- 
work of the dual color superconductor |^ |53 and in the 
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Gaussian Stochastic Model [Hill. The profile is rather 
well described by a Gaussian-like parameterization 



fgip) = Ng eyip[-lii2{p/pg 



(10) 



with a half width pg 0.3 fm and a steepness parameter 
ranging from n — 2.3 to n — 3.2. In tab.^we give two of 
the parameter sets used in 20], together with the given 
key quantities of the gg-string. 

All basic quantities of the Gornell potential and the 
profile of the gg-string agree either with lattice results 
or with results obtained in heavy meson spectroscopy, 
except for the Goulomb constant a which is somewhat 
small in our results. 



III. LONG RANGE BEHAVIOR OF gg-STRINGS 

In 01 we performed a detailed analysis of gg-strings 
with the main focus on the bulk properties of the strings 
such as the string tension r and the width of the strings 
Pg . The profile of the energy density, which we determine 
numerically by solving the Eqs. was well described 
by a generalized Gaussian form as in Eq. H10|l . On a lin- 
ear scale the deviations of the numerical results to the 
Gaussian fit were hardly seen. However, this parameteri- 
zation does not reproduce the long range behavior of the 
string fields far away from the string axis. Instead both 
the electric fields C and the confinement field a follow 
an exponential as can be seen in fig. El It should be noted 
that the exponential tail of the string fields does not in- 
fluence the results obtained in |23 • This exponential be- 
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FIG. 3: The logarithmic profile of the string. Shown is the 
profile of the confinement field x ~ ""vac — cr (solid line, 
cf. Eq. pi^ l and of the 3-component of the electric field D3. 
The slope of the D^-field is much steeper and is equal to that 
of the third power of the confinement field (dotted line) . 

havior can be explained quite naturally, if one assumes, 
that the electric fields die out on a much shorter char- 
acteristic length scale than the confinement field. With 
this assumption the source term in the right hand side of 



Eq. Ij3b|l vanishes. For small deviations 



X ■ 



(11) 



of the confinement field from its vacuum value we can 
make a Taylor expansion of U{<7) around dvac leading to 

U{a) « 5TOg(cr - CTvac)^- 

For strings with very large gg-separations R we can 
recast Eq. I|3b(l into 



ld_ 

pdp 



dx 



- m X = 



(12) 



where p is the coordinate transverse to the string axis and 
where we have used cylindrical symmetry. The regular 
solution of this equation is 



X{p) = XoKoimgp) w xo 



-nigp 



(13) 



where Kq (x) is the modified Bessel function of the second 
kind and the second relation holds for large p and xo is 
some constant. The parameter rUg is therefore directly 
connected to the screening mass mo- = rUg of the con- 
finement field or to its screening length ^ — l/rug. The 
electric field is = K{cr)E°' . In the chosen parameter- 
ization the function k((t) itself and its first two deriva- 
tives at a — (Tvac are proportional to the numerically 
small number Kvac • For small deviations x = o'vac — <J we 
expect 



D^ip) 



■niD 



«o(x/xo)^ + 0(kv 
«o ix/xof E 
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(14a) 

(14b) 
(14c) 



where kq is a proportionality constant and E is the elec- 
tric field which varies only slowly with p. The screen- 
ing mass of the electric displacement is therefore mu — 
3toct — 3mg and the field is screened on a characteristic 
length scale A — ^£,. This justifies a posteriori to ne- 
glect in Eq. (|12|l the source term on the right-hand side 
of Eq. (|3bl). We have added in fig. |3 the third power 
X'^ of X (dotted line) which shows the same slope as 
£)(a=3) (dashed line). To check this result numerically 
we have fitted the above analytical solutions in Eqs. I|13|) 
and 1)1 4bl) to both the confinement field x and the electric 
displacement Z?" obtained in our numerical calculations 
leaving the parameters mo- and rriD as fit parameters. In 
fig. 21 we show the results for different model parameters 
nig and for different string lengths R. The dependence 
of mo and m^i on rUg for fixed R is obvious. For fixed 
nig the fitted values approach roughly the expected val- 
ues rua- = nig and niD = 3m g, respectively, for growing 
string lengths R. Actually the fitted values undershoot 
the theoretical values consistently, which might be due to 
the final extent of the numerical box and the numerical 
boundary conditions x = on the box boundary. 
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No. 


fii/-* [McV] 


mg [McV] 


CTvac [fm~^l 


gs 




^ 1— 1 




Pa [fm] 


n 


I 


260 


1000 


1.29 


2.0 




980 


0.18 


0.33 


2.3 


1— 1 


240 


1500 


1.13 


1.8 


10-* 


980 


0.12 


0.34 


3.1 



TABLE II: In the first 5 columns we show the CDM parameter sets used in the description of gg-strings and gqg-baryons, in 
the last 4 columns we list the resulting values for the string tension r, the Coulomb-parameter a and the shape parameters of 
the profile, i.e. the width pg and the steepness parameter n as explained in the text. 
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FIG. 4: The fitted values of rricr (upper panel) and mo (lower 
panel) for different model parameters mg and for different 
string lengths R. 



IV. STRING INTERACTIONS 

In the CDM strings or flux tubes develop between two 
oppositely charged quarks, or more generally, between a 
number of particles with vanishing total net color charge. 
In the case of four particles, there must be two quarks 
and two anti-quarks to ensure color neutrality. For the 
discussion of string interactions the network of flux tubes 
does depend strongly on the configuration in coordinate 
space but also on the configuration in color space. In 
the latter there are two distinct possibilities. First, all 
quarks can be of the same color |c) and the correspond- 
ing anti-quarks have the anti-color |c), and thus the color 
content is (cccc). Second, there is one string with color 
content cc and another one with c'c' with c ^ c', so that 
the color content is {cc'cc'). This latter configuration is 
not possible in SU(2), where the two members of the fun- 
damental doublet are simultaneously the anti-particles to 
each other. In this work we concentrate on the former 
color configuration to compare our results to those ob- 
tained in SU(2) lattice calculations [sj Isol [60ll and in 
the model of the dual color superconductor |61| . 
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Both the confinement field a and the electric displace- 
ment D"^ follow an exponential far away from the string 
axis. They vanish on characteristic screening lengths 
with A = {3m g)^^ < £, = n^g^- This is somewhat sim- 
ilar to a dual type I color superconductor, where A and 
^ describe the penetration depth of the color field and 
the correlation length of the Higgs field. In ji^j and [s^ 
a lattice analysis of the string fields in Abelian projec- 
tion was made and a similar result A < ^ was obtained. 
However, it should be noted, that tjid — inig depends 
on the specific behavior of K{a) ai a — (Tvac and might be 
modeled differently. In contrast TOo- — mg follows from 
the e.o.m. as long as the ZJ^-field can be neglected in 
Eq. KM . 



FIG. 5: Upper panel; Four particles placed on the corners of 
a rectangle with length R and width D (left). Two parallel 
strings of length R and distance D (A). Lower panel; Two 
anti-parallel strings with R > D {B) and R < D {B'). Solid 
lines stand for the ground state of the two possible flux tubes, 
dashed lines for a possible excited state. 

In coordinate space the orientations of the flux tubes 
depend on the actual distributions of colors and anti- 
colors. We first discuss the simple case of four particles 
with colors cccc placed on the corners of a rectangle with 
length R and width D as shown in the upper left part of 
fig. El If all quarks q are lying on the left side and the 
anti-quarks q on the right side the ground state of flux 
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tubes will be like in configuration A in the same figure 
independent of R and D. In this configuration we have 
two strings of length R with the electric flux pointing in 
the same direction. In varying D we can examine the in- 
teraction energy of two strings of given length R and thus 
obtain the static string potential for this particular distri- 
bution. In (23l| the potential between two gg-strings was 
deduced with a dynamical setup, albeit with the oversim- 
plifying assumption, that the electric fields of two strings 
add coherently, thus neglecting the intrinsic four-particle 
interactions. 

If in contrast the positions of a quark and an anti-quark 
are interchanged, there exist two different possibilities. 
As the energy of a single gg-string is a monotonically 
rising function of separation R 'H?!, the lowest 

four-particle energy is obtained by minimizing the total 
string length. Consequently ii R > D, the flux tubes 
point upside down as in configuration B in fig.[Sl whereas, 
ii R < D, the flux tubes point from left to right as in 
configuration B' . In the CDM the ground state of the 
flux tubes is found within the model itself by solving the 
e.o.m. (|3J), without further input. If one increases D with 
fixed R, the strings will flip from configuration B to B' 
when D = R, and one can study the string flip interaction 
energy. In addition, a dynamical string flip can lead to 
the dissociation of a J/ip and thus to J/tp suppression in 
relativistic heavy ion collisions |31t] . 

We define the interaction between two strings as the 
difference between the 4-particle energy i?4 and the sum 
of the energies Es-^ + of two independent strings si 
and S2 with minimal energy. In the case of four parti- 
cles (7i(?2gi52 we calculate the sum of the energies Eq^q^ 
and £"^2 92 of the strings (gigi) and (9292) separately and 
compare it to the sum of the energies Eq-^q^ and Eq^q-^ of 
the competing pairs ((7192) and (92 In the absence of 
any interactions, the lower value defines the ground state 
energy of the strings si and S2 (solid lines in fig. El, the 
higher an excited state (dashed lines in fig. |SJ| . As the 
energy of a gg-string rises monotonously, we do not have 
to solve the equation of motion for both configurations 
but only for that with the minimal total string length. 
The interaction potential V4 is then given by 



V4 — i?4 — {Esi + Es^) , 



E., 



E, 



1112 



-^92 92 

^92 91 



(15a) 
(15b) 



equal. In that case we calculate the interaction energy 
density as 

Vi{r) = £4(r) - ^[(£gigi(f) +£g292(^) 

+ (£919-2 (^+ £929-1 (^)] , (17) 

i.e. we compare the four-particle density with an inco- 
herent superposition of the 2-string configurations B and 
B' . Of course this matters only in the description of the 
energy density and not of the potential energy V4. 

Given the exponential behavior of the strings far away 
from the string axis (cf. Eqs. (|13|l and one might 

estimate the interaction between two separated strings in 
the following way: Assume both strings point along the 
X-axis with their axes shifted by ±D/2 away from the 
string axis. For sufficiently large D, the approximations 
made in section IIIII are valid and the equation for the 
confinement field is linearized as in Eq. H12|l . The con- 
finement field of the two string configurations at |y| ^ Z? 
then simply is a linear superposition of two single strings, 
i.e. 



xiy) 

Xl/2 



xi{v) + X2{y) 
= t 



with 



(18a) 
(18b) 



with xo some constant. As the screening length A of the 
electric field is substantially smaller than that of the con- 
finement field, A < ^, the electric fields do not contribute 
much to the total energy density e, which consequently is 
dominated by the scalar potential U{x) ~ ^^^X^- The 
distribution of the interaction energy W4 (see Eq. (Uni)) 
can be expressed with Eq. (|18|l as 



My) = ^rnl{x^{y)-xliy)-xl{y)) 



^ — rriaD 



S^*^ TTlgD 



nigD 



(19) 



with Wo some constant and where we have used \y\ <C D. 
The strings at a distance D therefore interact with a 
Yukawa potential, which will be verified numerically later 



To illustrate, how the two strings of the ground state 
infiucncc each other we can define the spatial distribution 
of the interaction energy W4 as: 



Electric field D 



Viir) = ei{f) - (e^i (f) + e,^ (f)) 



(16) 



where the energy densities e are given as the integrands of 
Eqs. Ij4b|l - (|4d|) . One comment about the special config- 
uration where the particles are located on the corners of 
a square is needed. If the two strings are of type B/B' in 
fig-El both configurations are degenerate for R — D, that 
is the sums of energies Eq-^q-^ +Eq2q2 and Eq-^q.-^ ^Eq^q-^ are 



We start our numerical analysis by discussing the elec- 
tric field lines of two 4-particle configurations with square 
symmetry in fig. [HI where the particles are located at the 
corners of a square [R = D = Ifm). In this and the 
following figures we have a planar particle configuration, 
although the calculations are done in three dimension. 
The figures show a cut along the plane z = of the par- 
ticles. If not stated otherwise, the calculations for the 
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next figures were done using parameter set PSl from ta- 
ble^ In the upper panel, the strings are paraUel to each 
other (type A of fig. El and in the lower panel, the strings 
arc anti-parallel to each other (type B/B' in fig. |SJ). In 
both figures the field lines are shown both for the case 
K(r) — const = 1, i.e. for the color fields obeying the ordi- 
nary electromagnetic Maxwell equations (dashed lines), 
and for ^[(^(r)] as calculated from Eqs. l|3Jl (solid lines). 
All field lines are chosen to start on a circle around each 
quark (filled dots) with equal angular separation and con- 
sequently end on the anti-quarks (open dots). For both 
orientations the field lines in the CDM calculations are 
compressed into well defined flux tubes stretching from 
quarks to anti-quarks, as opposed to the k, — 1 case, 
where the field lines extend in all directions. The flux 
tubes of the anti-parallel configuration (lower panel) split 
symmetrically into two parts connecting each quark with 
both anti-quarks. Therefore we have a superposition of 
the two 2-string configurations B and B' of fig. [S] Due 
to the symmetry of the configuration, there is no electric 
field at all in the center at r = 0. 



B. Energy distribution 

To analyze the validity of our assumptions in Eqs. I|18(l 
and (|19|l we show the relative difference Ax = [x ~ (xi + 
X2)]lx (see Eq. (^HJ) for a configuration of two 1 fm 
long parallel strings in fig. [7| As expected, the difference 
approaches for increasing Z?, i.e. the scalar field in the 
center between two strings is a linear superposition of 
two isolated strings for sufficiently large D. Note, that 
Ax becomes concentrated at \y\ ^ D for large D, i.e. 
the fields of the strings cancel each other at the centers 
of each string. 

We expect that two strings of given length, that are 
asymptotically far apart from each other {D ^ i?), 
do only weakly interact. If they approach gradually 
[D K, i?), parts of the flux tubes overlap and both the 
scalar confinement field as well as the electric fields are 
distorted from their asymptotic shapes. We show the en- 
ergy densities e{r) of these distorted fields in fig.|Hl We 
choose a string length R = \ fm. In the upper panel, the 
two strings are parallel to each other (type A) and in the 
lower panel they are anti-parallel aligned (type B/B'). 
From left to right the string distance D decreases from 
D = 1.5 fm to £> = 0.5 fm. The symbols for the quarks 
are the white dots, those for the anti-quarks the black and 
white dots. The contour lines are equidistant in energy 
density in the range e = (0.4GeV/fm^ . . . 2.4GeV/fm^) 
in steps of Ae = 0.4GeV/fm^ 

For large separations {D = 1.5 fm > R) there are two 
nearly unperturbed strings independent of the orienta- 
tion of the strings. In this case the string configurations 
A and B' are indistinguishable. Only if the flux tubes 
get in contact with each other {D = 1 f m = R), the 
two orientations behave differently. For the parallel case, 
one finds only a slight attraction between the flux tubes 
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FIG. 6: The electric displacement D for k = k(o-) (solid lines) 
and for k = 1 (dashed lines). The strings have a length K = 
1 fm and their distance is Z) = 1 fm. The upper and the 
lower plot show a parallel and an anti-parallel orientation, 
respectively. The confinement mechanism of the CDM pushes 
the field lines into well defined fiux tubes. For the anti-parallel 
case the electric flux is distributed symmetrical into types B 
and B' . 



showing up in the distortion towards the center of the 
two strings. Also the energy density in the center of each 
string is lowered a little bit. 

However, for the case of anti-parallel strings, the flux 
splits up in the two directions and is the same superpo- 
sition of types B and B' as already seen in fig. We 
note, that this reorientation of the flux tubes is found 
just by solving the equations of motion (|3Jl. No exter- 
nal input such as the orientation of a Dirac-string as in 
the dual color superconductor model i6l| is needed. The 
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FIG. 7: The difference Ax = [x - (Xi + X2)]/x (cf. Eq. JUJ) 
for a type A configuration of two parallel strings of length 
R — 1 fm. The plot shows a cut along the y-axis transverse 
to the string axes through the center of the configuration. 



value K = 0.3 is characterizing the surface of the flux 
tube. The figure displays two different views of the same 
configuration, the first a perspective view of the system 
and the second a projection along the z-axis. Quarks and 
anti-quarks are marked as black and white spheres, re- 
spectively. The four lines are the shortest links between 
each quark/anti-quark pair. In this highly symmetrical 
configuration again the flux of each quark is split up into 
two equal parts pointing towards both anti-quarks, re- 
sulting in four equal flux tubes. The centers of the flux 
tubes are bent slightly towards the center of the configu- 
ration, as already seen before in the planar configuration. 
Moving one qq-string out of this symmetric configuration 
would cause two of the four flux tubes to break in favor 
of the other two. The string flip is therefore seen in the 
3-dimensional configurations as well. 



C. Interaction potentials 



CDM therefore incorporates the feature of string flip as 
used in the string-flip model ^62, ,65,] . For the situation of 
anti-parallel strings, the string flip is not a discontinuous 
process but a smooth transition. For large separations 
{D — 1.5 fm), there is already a small part of the energy 
flux that stretches to the transverse direction. The same 
is true for D = 0.5 fm with interchanged roles of types B 
and B' . The 4-particle flux tube network is therefore a 
superposition of types B and B' with continuously vary- 
ing relative strength. 

If the two strings approach even further {D = 0.5 fm < 
R) the two parallel flux tubes melt into an extended flux 
tube and form a single bag (fig. El upper right panel). We 
observe the contact of the two pairs of likewise charged 
particles which will lead to a repulsive Coulomb interac- 
tion. For the anti-parallel situation however, the string- 
flip from type B' to B has nearly completed and the flux 
tubes point upside down (flg.|Hl lower right panel). The 
flux tube pattern now shows two flux tubes of length 
D < R. Because of this, the deflnition of the strings is 
somewhat ambiguous. We start in varying the distance 
D between two i? = 1 fm long gg-strings and flnally end 
up with two gg-strings of length D. Note that in this 
whole section we denote with the string length the par- 
ticle distance that we keep fixed, if not stated otherwise 
explicitly. 

Of course, four particles do not have to be in a plane 
necessarily. As our numerical realization does all calcu- 
lations in three dimensions we can easily describe flux 
tubes stemming from any arbitrary 3-dimensional parti- 
cle configuration. As one example we show the flux tubes 
of four particles placed on the corners of a tetrahedron 
with pairwise distance i? = 2 fm. We choose this large 
particle separation to clearly show the resulting flux tube 
network. In the 3-dimensional illustration in flg. |51 the 
equipotential surface of the dielectric function k = 0.3 
is shown. In the center of a gg-string k « 0.8 for both 
parameter sets given in tab. ^ (see j^). Therefore the 



The energy of a gg-system scales with the particle 
distance R according to the Cornell-potential given in 
Eq. (O. In the previous section we have shown the en- 
ergy distribution of two such interacting strings as well 
as the difference of the 4-particle state to the incoher- 
ent superposition of two equivalent qq-strings. From the 
integral of the energy density we get the total energy 
E4 in Eq. Q and from the corresponding difference of 
energies in Eq. (|15|l we extract the interaction potential 
V4. For the following calculation we keep the individual 
string length R and the relative orientation flxed, and 
vary only the distance D between the string centers. 

The total energy E4 of two i? = 1 fm long flux tubes is 
shown in flg. 1101 The orientation is parallel in the upper 
panel and anti-parallel in the lower panel. The total en- 
ergy (solid line) is separated in the different energy parts 
according to Eq. (0J. Here we have subtracted from the 
total energy the energy Eao of two inflnitely separated 
qq-strings. In the parallel case this is identical to the 
potential V4, as there is no string flip. For clarity we 
have separated the curves for the three energy fractions 
by equal offsets of 0.1 GeV (upper panel) and 0.5 GeV 
(lower panel), respectively. 

The total energy saturates for distances D > 1.5 fm, 
i.e. the strings cease to interact. If they approach each 
other, we flnd for the parallel strings a stable distance 
at D 0.4 fm. For smaller distances the two like- 
wise charged particles experience the Coulomb repul- 
sion, which is seen only in the electric part of the energy 
(dashed line). The volume part of the energy (dashed- 
dotted line) exhibits a small maximum when the two flux 
tubes get in contact. When they approach each other, the 
tails of the strings overlap and the two-string configura- 
tion takes on a larger effective volume than two isolated 
strings. As a consequence, the volume energy rises while 
the electric energy drops down. 

For the anti-parallel strings the energy behaves similar 
for D > 1 fm. However, for smaller distances the string 
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FIG. 8; Energy density of two _R = 1 fm long gq-strings for different distances D. The orientation of the strings is parallel 
(upper panel) and anti-parallel (lower panel). Contour lines are chosen at the values e = (0.4 GeV/fm'^ . . . 2.4 GeV/fm^) in 
steps of Ae = 0.4GeV/fm^. The darker the area the larger is e. Quarks are depicted as white dots, anti-quarks as black and 
white dots. Dark shaded areas are due to the strong Coulomb fields at the particle positions and the contouring is left out here 
for better visibility of the flux structure. 



flip takes place and the energy behaves nearly as two flux 
tubes that shrink with D. 

In fig. we present the interaction potential V4 for 
different string lengths R. The orientation is parallel 
in the upper panel and anti-parallel in the lower panel. 
The generic form does not change with R, but both the 
position of the minimum and the depth of the poten- 
tial change with R. For the parallel strings the min- 
imum of the potential is moving with increasing R to 
larger separations D. For string lengths R exceeding 
1 fm the location of the minimum basically stays con- 
stant. It reaches a stable point at 13 « 0.4 fm which 
is roughly the same size as the radius of the single flux 
tubes (see tab. Qualitatively this potential resembles 
the nucleon-nucleon interaction with a short/long range 
repulsion/attraction and a dip of the order of 100 MeV. 

The potential of the anti-parallel strings exhibits a 
characteristic kink at D = i?, which is due to the string 



flip. At this point the orientation of the strings flip from 
configuration B to B'. The results obtained in this work 
differ both qualitatively and quantitatively from that in 
[2^. In the older work it was assumed, that the electric 
fields of the to strings might be added linearly, thus ne- 
glecting the 4-body interactions of the color charges. In 
addition, the parameters chosen there led to flux tubes 
with a radial width of « 1 fm, which is much larger 
than our value pg « 0.4 fm (cf. tab. EJ. Consequently, 
the string interaction in the work of Loh et.al. has a range 
of about 1 fm. Also there was no stable point in the po- 
tential as is seen in our work in fig. 1111 

In fig. 1121 we have isolated the potential minimum for 
different string lengths R. In the parallel case (open 
triangles) the minimum scales linearly with the string 
length. This can be understood, as for long parallel 
strings the profile does not change along the flux tube 
axis. Thus the energy gain per string length should be 
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FIG. 9: (color online). Four particles on the corners of a 
tetrahedron. All pairwise distances are equal to i? = 2 fm. 
Shown is the equipotential surface n{f) — 0.3. Quarks and 
anti-quarks are marked as dark and light spheres, respectively. 
The lines are the direct links between each quark/anti-quark 
pair and are given to show the distortion of the interacting 
flux tubes. The upper and the lower plot give different views 
of the same configuration. 



constant when the two strings get in contact. 



For the anti-parallel configuration aX D — R but also 
for the tetrahedron configuration, when all pairwise par- 
ticle distances are the same, we can parameterize the 
4-particle potential in the spirit of the 2-particle Cornell 
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FIG. 10: Energy £4 of two parallel gg-strings (top) and two 
anti-parallel gg-strings (bottom) with length R = 1 fm. The 
total energy (solid line) is split into the different energy parts. 
All energies are reduced by their energies Eoo of two strings 
asymptotically far apart. The energy parts are separated by 
constant offsets of 0.1 GeV (upper panel) and 0.5 GeV (lower 
panel) for better visibility. 



potential VdR) ®: 



R V2R 



CFa + ATiR-2Vc{R) 



-(2-V2] ^ + (4r4 - 2t) R 



R 
4 2 

'r~~r 



(4r4 - 2r) R 



(20a) 

CFa + ATiR-iVciR) 

(20b) 



Here the first term in each bracket is due to the four 
attractive gg-pairs and the second to the two pairs qq 
and qq. The constant terms cancel each other exactly. 
The Coulomb interaction is only partially reduced in 
the square configuration and completely canceled in the 
tetrahedron geometry. 

In the above parameterization r4 is an effective 4- 
particle string tension, whereas r denotes the standard 
gg-string tension from Eq. @. Naively one might esti- 
mate r4 in the following way. As we have seen in fig. |S1 
(lower middle panel) and in fig. |51 the fiux tubes of the 
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FIG. 11: Interaction potential 14 of two parallel gg-strings 
(top) and two gg-strings (bottom) with different lengths R. 
The kink in the anti-parallel configuration is due to the string 
flip aX D = R. 



anti-parallel square configuration and of the tetrahedron 
configuration do not overlap for string lengths R> 1 fm. 
The electric flux from each quark is therefore split into 
two equal flux tubes pointing to the two anti-quarks. 
It was shown in [23|, that in the CDM the string ten- 
sion scales with gs^/Cp, i.e. with the charge of the parti- 
cle. For the two symmetric configurations the 4-particle 
string tension therefore should be half of the qq-string 
tension = t/2. In this case the linear confinement 
terms in Eq. 1)2(1 al) and Eq. (|20b|) vanish as well. 

In this simplified picture, the potential Vu of the pla- 
nar anti-parallel configuration should scale Coulomb-like 
and that for the tetrahedron should vanish very rapidly, 
compared to the former one. The result of the CDM cal- 
culations is shown in fig. E| For both the anti-parallel 
configuration (solid squares) and the tetrahedron (solid 
triangles), the potential Vi^uiin scales linearly with the 
string length for R > 1 fm. Neither a Coulomb-like nor 
a rapid cancellation is observed. Thus the 4-particle po- 
tential is not a trivial combination of (/^-potentials in the 
sense of the generalized Cornell potentials as given in 
Eqs. (EOJ. 

The string-string potential was also analyzed in SU(2) 
lattice theory for anti-parallel configurations in [s^l . The 



0.5 1 1.5 2 

R [fm] 

FIG. 12: The minimum of the potential 14(1)) for differ- 
ent string lengths R. The magnitude of the potential Vi.min 
rises linearly in the parallel case (open triangles) as expected 
for long homogeneous strings. In the anti-parallel case (solid 
squares) and the tetrahedron case (solid dots) the attraction 
gets weaker and shows an unexpected linear increasing be- 
havior. 



qualitative behavior of the potential is the same as in our 
model, although the absolute values of the potential V4 
are consistently smaller than ours. However, it should 
be noted that the parameters of our model used in this 
work are not adjusted to the 4-quark problem, but are 
fixed to the gq-properties only. In a very recent SU(3) 
lattice calculation of the 4-particle system js^ a multi- 
Y flux tube picture was proposed, such that the total 
string length is minimized. This is similar to the Y-like 
flux tube picture of the 3-quark system. In our model 
this would be possible with another color content like 
e.g. rfgg, which is devoted to future work. 

Next we turn to the long range behavior of the poten- 
tial ¥4(0). From Eqs. p9|l we expect a Yukawa potential 
for sufficiently far separated strings. In fig. 1131 we show 
the negative of the potential V4 for different string lengths 
i? in a parallel orientation with parameter set PSl. The 
dotted lines are the best fits of a Yukawa type poten- 
tial Vyuk{D) = Vq exp {—nisD) /{rUsD) as expected from 
Eq. itT^ to the CDM potentials, with rus and Vq being 
fit parameters. For D > R the potential follows nicely 
the Yukawa behavior. From Eqs. H19|) we expect also, 
that the screening mass TOs is given by the curvature nig 
of the scalar potential U at a = (Tvac- To test this we 
show the dependence of on mg in fig. 1141 From top 
to bottom, the relative orientation of the strings is par- 
allel, anti-parallel and transverse (tetrahedron like). We 
have extracted the screening mass from the fit for dif- 
ferent string lengths R and for the two parameter sets 
PSl and PS2 from tab. ^1 The parameter nig takes on 
the value rrig = 1000 MeV (PSl) and rrig = 1500 MeV 
(PS2), respectively. The error bars shown in the plot 
result from variations of the fit interval of D used in 
the fit. We note, that 14(1)) is difficult to extract nu- 
merically, as it is an exponentially small difference be- 
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FIG. 13: Logarithmic plot of the potential V4,{D) for different 
string lengths R. The dotted lines are fits of a Yukawa type 
potential to V4{D). The model parameter correspond to PSl 
and the strings are parallel to each other. 



tween to large energies (see Eq. (|15|) '). It should be 
noted also, that it is numerically more difficult to cal- 
culate the potential V4 for large string lengths R due 
to the limited numerical box size but also for parame- 
ter set PS2 due to the strongly pronounced maximum 
of the potential U (see fig. EJ. Therefore the error bars 
are larger for longer strings but also for parameter set 
PS2. However, within the errors a good agreement of 
the screening mass to iris = rUg is seen. This behavior 
is almost not dependent on the string length R but also 
not on the relative orientation. In the dual Ginzburg- 
Landau model a Yukawa-type potential was found with 
a screening mass = 1430 MeV IgJ, which is consis- 
tent with our results. Following our above interpretation 
of rUg , one might compare these numbers to the glueball 
mass which has been calculated on the lattice between 
1500 MeV <mg< 1700 MeV [Hil. Another possibil- 
ity is to compare rUg with the mass WoS ~ 1200 MeV of 
the off-diagonal gluons given on the lattice as well [isf . 
The detailed verification of the Yukawa potential between 
two strings, as proposed in our description via the CDM, 
should be a task for future lattice calculations. 



V. MULTI-QUARK SYSTEMS 

In this section we present CDM results of overall color- 
less multi-quark systems, with the particle number being 
larger than four. Such states might in principle exist and 
a number of possibilities like the pentaquark j68i . <6^ , the 
H-dibaryon jTfl and strangelets |71h7^ , were eagerly dis- 
cussed in the literature. Moreover we want to look for a 
possible transition if the quark number density becomes 
large as e.g. in the interior of dense neutron stars or in 
the highly compressed phase in relativistic ion collisions. 
We study unordered ensembles of quarks and anti-quarks 
with a varying number of particles. The particles are 
placed in a given volume which is sufficiently smaller than 
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FIG. 14: The screening mass of the interaction potential be- 
tween two strings for PSl (solid error bars) and PS2 (dashed 
error bars). From top to bottom the strings are oriented 
parallel, anti-parallel and transverse (tetrahedron like) to 
each other. Within the uncertainties rUs is equal to nig = 
1000 MeV (PSl) and nig = 1500 MeV (PS2). 



the numerical box in order to reduce boundary effects. 

The color of the particles are chosen according to 
two different schemes. In the first we first pick a color 
c S {r, (7, 6}. This color is assigned to a quark q and si- 
multaneously the corresponding anti-color c is assigned 
to an anti-quark q. This gg-pair is thrown randomly 
into the volume. We repeat this procedure until a given 
number Nm of qq-statea is reached. In this way we get 
a system with vanishing net baryon number, although 
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FIG. 15; (color online). Flux tube structure of systems with baryon number TV;, > (upper panel) and with vanishing baryon 
number (lower panel) for different particle densities. We show the equipotential surface of the dielectric function ki, — 0.4 
characterizing the surface of the flux tubes. Color codes for the particles are red/green/blue and cyan/magenta/yellow for the 
quarks and anti-quarks, respectively. The particle numbers for the upper row from left to right are A'^ = (63, 126, 255) and for 
the lower row A'' = (64, 128, 256) corresponding to particle densities n = (0.5, 1.0, 2.0) fm~^. 



baryonic and anti-baryonic clusters might be formed. In 
the second scheme, we assign the three colors r,g,b to 
three quarks only and throw them into the volume un- 
til a given number N}, of baryonic qqq states is reached. 
In this scheme the baryon number is Nf, > 0. Finally we 
vary the total number of quarks N = 2Nm and N — 3Nf,, 
respectively, to study the behavior of the quark system 
for different particle densities n — N/V. For each particle 
density n we choose many different spatial configurations 
to calculate average quantities as the energy per particle 
and the bag volume per particle. For this rather quali- 
tative analysis we restrict ourselves to the parameter set 
PSl. For percolation studies, such as the quark density 
dependence of the mean bag size, or the formation of a 
single super cluster, we would need more statistics. This 
analysis is devoted to future work. 

In fig. El we show the resulting flux tube structures 
for the baryonic (upper panel) and the mesonic case 
(lower panel), respectively. In the latter one one can 
find all possible color neutral subsystems, namely qq- 
states, qqq-states and qqq as well as subsystems with 
larger numbers of quarks and anti-quarks. For these fig- 
ures the particles were thrown into a cubic box with a 
size V — (5fm)^. To visualize the flux tubes in three- 
dimensional space, we show the position of the particles 
as small spheres and the equipotential surface of the di- 
electric function at a value k = Kf, = 0.4 characterizing 
the surface of the flux tubes or bags. The number of 
quarks in the baryonic case (upper panel) in fig. 1151 from 



left to right is = 63, 126, 255 corresponding to parti- 
cle densities n = (0.5, 1.0, 2.0) fm~^ or baryon densities 
rib = n/3 = (1.0,2.0,4.0) x no with no = 0.17 fm"^ be- 
ing nuclear matter density. The number of quarks and 
anti-quarks in the mesonic case (lower panel) from left 
to right is iV = 64, 128, 256 leading approximately to the 
same particle densities as in the baryonic case. 

The formation of well defined bags is clearly seen. For 
small particle numbers the system is dominated by long, 
nearly linearly shaped flux tubes. In the baryonic en- 
semble the smallest clusters are ggg-states, whereas in 
the mesonic ensemble the smallest clusters are in prin- 
ciple gg-states, ggg-states and ggg-states. For denser 
systems the average number of particles per bag grows 
but still distinct isolated bags are formed. Even for the 
n — 1.0 fm~^ systems (the two centered figures) the flux 
tube structure is maintained. The individual and still 
separated forms are very complex objects and can be in- 
terpreted as multi-quark excitations of hadronic particles 
much akin to the old bootstrap picture of higher lying 
resonance states {7'k i75j . The non-perturbative vacuum 
is replaced by a spaghetti like perturbative vacuum and 
both vacua still balance each other. In the systems with 
n = 2.0 fm^^ (the two figures at the right) nearly all 
particles are gathered in one single but highly deformed 
bag. The transition from a small particle density with 
distinct bags and a small number of particles per bag to 
a large particle density with only one super-cluster and 
a large number of particles per bag is similar to a per- 
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eolation transition which was proposed to occur for the 
quark hadron transition |76l - l7a | . 

In the following we analyze this interplay between 
the perturbative and the non-perturbative vacuum more 
quantitatively. To reduce numerical boundary effects, we 
reduce the volume, where the particles are thrown in, to 
a spherical region with volume V = |7rr^ and radius 
r = 2.5 fm. In fig. 1161 we show the energy per particle 
as a function of particle density in a double logarithmic 
plot. The upper and the lower plot show a baryonic and 
a mesonic system, respectively. The solid squares are 
the values for the total energy per particle averaged over 
many configurations for each density. The error bars de- 
note the statistical standard deviation for the ensembles 
of configurations. Clearly, the statistical fluctuations of 
the energy is largest for the smallest densities. For com- 
parison, we show the electric part of the energy (solid 
triangles) and the equivalent energy per particle Eq/N 
for the pure Maxwell case, i.e. for k = 1 everywhere 
(open dots). The fluctuations for the energy per par- 




FIG. 17: (color online). 654 quarks randomly distributed in 
a sphere with radius r — 2.5 fm and volume V. The non- 
perturbative phase (cr — (Jvac) is pushed out of the sphere 
completely and Vport ~ V. The bag surface is defined by 

K — Kb ~ 0.4. 




n [1/fnT'l 

FIG. 16: The energy per particle E/N according to Eq. ||1J 
as a function of particle density n for a baryonic (top) and a 
mesonic (bottom) system. Solid squares symbolize the aver- 
age values at each density with the statistical standard devi- 
ations denoted by the error bars. Triangles are the average 
values of the electric energy only and the open dots the av- 
erage energy per particle for the Maxwell case (k = 1). The 
dashed line is a fit to the average energy per particle. 

tide are largest for small densities n, as the particles 



get more and more homogeneously distributed in space 
with growing density. Both the average total energy and 
the electric energy follow roughly a power law, whereas 
the free energy Eq/N stays constant. We can estimate 
the scaling of the total energy as follows. The energy 
of a cluster scales linearly with the size L of the clus- 
ter. For a given particle configuration, the clusters form 
themselves by minimizing the total energy. Each colored 
quark builds up a flux tube to the nearest oppositely 
colored anti-quark or sub-cluster. Thus we can estimate 
the cluster size by the average particle distance of the 
system, i.e. L = n^^/"^. The total energy per particle 
is therefore Etot/N oc n~^^^. We have fitted the to- 
tal energy per particle to the ansatz Etot/N — cn^ in 
the range n < lfm~^. The result for the fit parame- 
ter was (3 — —0.35 and (3 = —0.36 for the baryonic and 
the mesonic system, respectively, which is close to the ex- 
pected result (3 = —1/3. The fit to the low density region 
overshoots slightly the high density results (n > 3fm~'^), 
indicating that the assumption of isolated long flux tubes 
is not valid anymore. Instead, for very large densities, the 
whole volume is filled with particles homogeneously and 
the non-perturbative phase is pushed out of the volume 
V . The dielectric function is therefore n < 1 everywhere 
inside V and the electric energy should be the same as 
for the equivalent free case. Indeed, the electric part of 
the energy (triangles) approaches the free energy (open 
dots) for densities n > 5fm^'^. 

Another interesting quantity is the specific bag vol- 
ume per particle Wport = ypcrt/N. For small n one finds 
dominantly bags with two or three particles. The size of 
the bags per particle decrease with the average particle 
distance, i.e. with increasing n, until with still further 
increasing density the bags start to overlap and melt to- 
gether. The decrease of the bag volume per particle slows 
down or might even turn over to an increase with increas- 
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ing n in the same way as we have seen in the bump of 
the volume energy for the melting gg-system in fig. IIUI 
At a critical density ric, when the melting of the bags 
is completed and the non-perturbative phase is pushed 
out of the volume V completely, the total perturbative 
volume can not increase anymore and Vport = const = V 
as shown in fig. 1171 With further increasing density f port 
decreases as n^^. This critical density ric marks the tran- 
sition to a system of deconfined quarks and anti-quarks. 



10" 



> 10"' 



10"' 







baryonic 


K=0.6 — ^-! 


1 




if H*' 

















0.1 1 
n [1/fm^l 



10" 



> 10"' 



10"' 







mesonic 


K=0.4 ' ■ ' 
K=0.6 ^--J.-! 






1 


4 {III'. 




I 















0.01 



0.1 



10 



n [l/fm-^l 



FIG. 18: The specific bag volume per particle for a baryonic 
(top) and mesonic system (bottom). Squares denote the av- 
erage specific bag volume for Kb — 0.4 and triangles that for 
Kb = 0.6. Error bars indicate the standard deviation for the 
set of different configurations at each density. At Uc ~ 3fm~^ 
the melting of the single bags is completed and the quarks are 
in the deconfined phase. 

We measure the total bag volume as that part in space, 
where k > Kb. To show the dependence of the bag volume 
on Kf, we choose the two values kj, = 0.4 and ki, = 0.6, 
where Kb = 0.4 (Kb — 0.6) is a surface somewhat more 
in the exterior (interior) of the bag. Therefore the bag 
volume measured in this way is larger for Kb — 0.4 than 
for Kb = 0.6. The result is shown in fig. ^| for the bary- 
onic (top) and the mesonic case (bottom) for Kb = 0.4 
(squares) and Kb = 0.6 (triangles). Again the error bars 
indicate the standard deviation of the specific volume 
for the different measured configurations at each density. 
The bag volume per particle shows exactly the antici- 
pated result. It decreases with increasing n for small den- 



sities. For 0.3 fm~^ < n < 3 fm~^ the specific bag volume 
Wpcrt develops a plateau for Kb = 0.4 and rises again for 
Kb = 0.6. We find the critical density ric, i-e. the point 
where Wpcrt decreases again, at around ric ~ 3.0 fm""^. 
This corresponds to a baryon density rib — ric/S ~ 6no, 
i.e. six times nuclear matter density or a meson density 
of approximately = ric/2 — 1.5 fm^'^. We note, that 
even at the highest density n — 10 fm^'^ the mean parti- 
cle distance is much larger than the intrinsic particle size 
To of the particles, so that the color charges do not cancel 
out each other by overlapping. 
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FIG. 19: The balance between the electric and the volume 
energy for a baryonic (top) and a mesonic system (bottom). 
The open squares and the solid dots denote the values of ev- 
ery configuration and their average value, respectively. The 
two energies are the same within 10% for the numerically rea- 
sonable range. 

It was one of the results of the work in j20[] , that the 
volume energy E^oi and the electric energy Ec\ of a qq- 
string nearly balance each other. This holds for larger 
particle densities as well until the system reaches the 
deconfinement transition. In fig. 1191 we show the ratio 
Evoi/ Eel as a function of the particle density for a bary- 
onic (top) and a mesonic system (bottom). It is constant 
within 10% (neglecting the first point with large uncer- 
tainties in the mesonic system) for densities n < 3.0 fm""^ 
but falls down for higher densities once the system is in 
the deconfined phase. The bag volume and therefore the 
volume energy i?voi are maximal but the electric self- 
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FIG. 20; Energy density as a function of particle density n. 
Baryonic (squares) and mesonic systems (triangles) are indis- 
tinguishable from each other. The energy density is well de- 
scribed by a power law (dashed line) and is a smooth function 
even across the critical density at ric = 3.0 fm""^. The verti- 
cal and horizontal lines mark the uncertainties in the critical 
particle density and critical energy density, respectively. 



energy associated with the particles increases with the 
particle density. 

Finally we show the total energy density e = i?tot /V as 
a function of the particle density in fig.|2n|for both bary- 
onic (squares) and mesonic (triangles) systems. It fol- 
lows nicely a polynomial behavior (dashed line) over the 
whole density range. Remarkably, it is indistinguishable 
for (static) mesons and baryons. At the critical density 
He = 3.0fm~'^ it has a value Sc = 1.2GeV/fm'^. Due to 
the uncertainty in the exact value of the critical density 
Uc we give a band for the critical energy density between 
l.OGeV/fm^ and 1.3GeV/fm^ This is in good agree- 
ment with the critical energy density found on the lattice 
[t^. where it was found to be between 0.5 GeV/fm'^ and 
l.OGeV/fm^ 

We note, that there are basically no differences in 
figs. El a-iid El between baryonic and mesonic sys- 
tems. Apparently all effects are driven by the particle 
density, no matter which kind of particles (quarks and 
anti-quarks or only quarks) are studied. 



VI. SUMMARY AND DISCUSSION 

The description of the interactions of quarks and glu- 
ons, both at high temperatures/densities and in vac- 
uum, from first principles is still an open task to solve. 
The chromodielectric model offers a tool to describe in 
a transparent way these interactions while including the 
confinement phenomenon dynamically within the same 
framework. 

Within the chromodielectric model particles group 
themselves into color neutral bags with the means of color 
electric flux tubes. In this work we studied the interac- 
tions of these flux tubes with each other. We have seen 



both on the level of the electric fields and on the level of 
the energy density, that the flux tubes attract each other. 
Depending on the relative orientation of the strings, the 
flux tubes either melt gradually together or change the 
direction of the electric flux via a characteristic string 
flip. This should also be seen in future lattice SU(3) cal- 
culations for many heavy quarks. 

The attraction is seen also in the interaction potential 
between the strings. The depth as well as the range of the 
interaction changes with the length of the strings. The 
interaction depth in the order of a few hundred MeV is 
relatively strong compared to a rather old analysis made 
in lattice SU(2) calculations but also compared to a sim- 
ilar analysis within the dual color superconductor model. 
In the case of the string flip situation the potential can 
not be described by a simple incoherent superposition 
of flux tubes, but shows a real multi-particle effect. The 
long range behavior of the potential can be well described 
by a Yukawa type potential. The exponential decay of the 
interaction is expressed by a screening mass between the 
flux tubes which is rather insensitive to the relative orien- 
tations of the strings and to their length. The screening 
mass scales with the mass parameter mg introduced in 
the scalar self-interaction of the confinement field. It is a 
future task to check the model continuously against up- 
coming results performed on the lattice. On the other 
hand the Yukawa interaction should be tested in detailed 
future four-quark lattice calculations. Only recently we 
have become aware of a lattice study about the interac- 
tion of four- and five-quark systems ,34 . .80.] . 

In the studies of the multi-particle systems we have ex- 
plored the structure of the corresponding flux tubes. The 
non-perturbative phase is pushed out of the system with 
increasing particle density. Even for large particle densi- 
ties compared to nuclear matter constituent quark den- 
sity, the systems show a heterogeneous structure rather 
than a homogeneous transition. This situation resem- 
bles a typical percolation transition from a hadron to a 
quark gas. At low particle densities the energy scaling 
with the particle density is characteristic for a system of 
strings whose energy scales with the size of the strings. 
By increasing the particle density, the flux tubes or quark 
bags start to overlap and to melt into each other. The 
transition region between isolated hadronic objects and a 
single large super-cluster, where the quarks behave as de- 
confined particles, is found to be between n = 1 fm~'^ and 
n = 3fm~'^. The critical energy density of the transition 
to the deconfined phase is given between e = l.OGeV/fm'^ 
and s = 1.3GeV/fm^. This transition shall be examined 
in the future more carefully in a more detailed perco- 
lation analysis. There one could measure for example 
the onset of the super-cluster, or the distribution of the 
bag size with respect to its volume or to the number of 
particles belonging to it, all as a function of the particle 
density. 
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